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ABSTRACT 

We present a study of internal stellar population gradients in early-type cluster galax- 
ies. Using the VLT VIMOS integral field unit, we observed 19 galaxies in the core 
of the Shapley Supercluster (z = 0.048). The radial trends in nine absorption lines 
(HSF to Fe5406) were measured to the effective radius for 14 galaxies, from which 
we derived the gradients in age, total metallicity and a-element over-abundance. We 
combine these with results from 11 galaxies studied in our previous VIMOS work 
(jRawle et al.ll2008b [). We observe a mean metallicity gradient of -0.13 ± 0.04 dex~^ 
and, in common with the findings of previous studies, galaxies with log cr ^2.1 have 
a sizeable intrinsic scatter in metallicity gradient. The mean log(age/Gyr) gradient is 
-0.02 ± 0.06 dex~^, although several galaxies have significant positive or negative age 
gradients. The mean gradient in a-element enhancement is -0.10 ± 0.04 dex~^. 

We find that stellar population gradients are primarily related to the central 
metallicity: early-type galaxies with super-solar centres have steep negative metal- 
licity gradients and positive age gradients; those with solar metallicity centres have 
negligible [Z/H] gradients and negative age gradients. There is a strong observed anti- 
correlation between the gradients in age and metallicity. While a part of this trend can 
be attributed to the correlation of measurement errors, we demonstrate that there is 
an underlying intrinsic relation. For the Shapley galaxies, B-R colour gradients pre- 
dicted from spectroscopic age and metallicity generally agree well with those measured 
directly from photometry. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: stellar content 



1 INTRODUCTION 

For almost four decades, observations have identified sys- 
tem atic radial tren ds in the stellar properties of galax- 
ies dSandagiEazS)- 

These internal gradients are a valu- 
able probe into the star formation histories and chemical 
evolution of galaxies. Traditionally, the focus has been on 
radial trends in (mostly optical) broadband colours (e.g. 
iPeletier et allllQQOl : pTamura et al.ir2000l : lMichardll2005l ). 

Early-type galaxies typically become steadily bluer 
towards their outer regions. These colour gradients 
have often been interpreted as variations in metallicity: 
iTamura fc Ohtal (|2003) estimated that VB-R ~ -0.1 corre- 
sponds to V[Z/H] ~ -0.^]. The relative shallowness of these 
derived metallicity gradien ts rule out th e simplest models of 
classic monolithic collapse (|Larsonlll974h . Within the frame- 



* Based on observations from tfie European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile (pro- 
grammes: 078.B-0539, 081.B-0539) 
f E-mail: t.d.rawle@durham.ac.uk 

^ All gradients are quoted per decade radius, with VX the short- 
hand for the change in X per decade radius 



work of the hierarchical ACDM cosmology, competing pro- 
cesses influence the gradients. For instance, a galaxy with 
a history of m erging would show diluted or erased metallic- 
ity gradients l|White|[l980h . but a central star burst could 
strengthen the negative metallicity gradient and create a 
positive age gradient (|Barnes fc H crnguist 1991'). 

The well-known degeneracy between age and metallicity 
in broadband optical colours has hindered attempts to dis- 
tinguish their respective contribution to the gradients. Using 
optical an d infrared colours to estimate metallicity and age 
gradients, iTamura &: Ohtal (|2004i) found that including the 
age effect weakens the mean V[ Z/H] (to -0 . 16 ± 0.09), while 
predicting small age gradients. IWu et al] (|2005| ) used opti- 
cal/IR colour gradients of 36 nearby early- types to derive 
a mean metallicity gradient of -0.25 ± 0.03 and a shallow 
mean Vlog(age) = 0.02 ± 0.04, but with a substantial dis- 
persion in age gradients of 0.25. 

Spectroscopic data can be used to discrimin ate more 
clear ly between the age and metallicity effects l|Worthevl 
Il994 ). To date, most spectroscopic studies of the radial 
trends of early-ty pe galaxies have used a single slit along the 
major axis (e.g. iGorgas. Efstathiou fc Aragon Salamancal 
Il990l : iMehlert et al.l l2003l : ISanchez-Blazguez et all l2007l : 
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iReda et all |2007| ). occasionally supplemented by a per- 
pe ndicular slit along the minor axis (e.g. the SO study 
in iNorris. Sharpies fc Kuntschnedl2006t l. Long-slit spectro- 
graphs have good optical throughput, but a significant frac- 
tion of the galaxy light goes undetected, and the major axis 
is assumed to be representative of the whole. 

An integral field unit (IFU) couples the diagnostic 
power of spectroscopic observations, to the spatial sampling 
enjoyed in imagin g. T he SAURON project, introduced in 
iBacon et al] (1200 ll ) and lde Zeeuw et al.l l|2002h . observed 72 
nearby galaxies {cz < 3000 km s~^). The excellent spatial 
coverage of the instrument permitted two-dimensional dy- 
namical/kinematic mapping of the targets. However, the 
trade-off for the spatial resolution is a limited spectral 
range covering only three useful absorption line indices (H/3, 
Fe5015, Mgb51 77), hindering the det ailed analysis of stel- 
lar populations. iKuntschner et al] (|200 6) showed that while 
the spatial maps of H/3 suggest that early-type galaxies typ- 
ically exhibit no radial trend in age, Fe5015 indicates that 
metallicity gradients are generally negative (mean VFe5015 
= - 0.014). 

iRawle et al] l|2008bl . R08) employed the Very Large 
Telescope (VLT) Visible MultiObject Spectrograph (VI- 
MOS) IFU to study galaxies in the local cluster Abell 
3389 {cz = 8100 km s-^). The VIMOS field-of-view is weU 
matched to the size of early-type galaxies at z = 0.02 - 0.05, 
and the instrument permits a large optical wavelength range 
covering many useful absorption line indices. Rather than 
creating two-dimensional maps of galaxies, radial trends of 
the spectroscopic features were measured out to the effec- 
tive radius. From these, gradients were derived in the stel- 
lar population parameters (age, metallicity and a-element 
abundance). Generally, galaxies displayed negative metal- 
licity gradients and flat age gradients. 

The present paper expands the sample of R08 with new 
observations of 19 early- type galaxies in the Shapley Super- 
cluster. The Abell 3389 targets have been re-analysed using 
the improved techniques described below. 

This paper adopts the cosmological parameters 
{nM,O.A,h) = (0.3,0.7,0.7). At the redshift of the Shapley 
Supercluster {z = 0.048), 1 arcsec corresponds to 0.95 kpc. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample and observations 

The sample consists of 19 early-type galaxies located in 
the three main constituent clusters (Abell 3556, Abell 3558, 
Abell 3562) of the Shapley Supercluster core {<cz> « 14660 
km s~^). Shapley is the riche st supercluster in the nearby 
universe (lRavchaudhurviri989l ) . and there is a wealth of an- 
cillary data for each galaxy. Photometric coverage includes 
ESQ 2.2m B- and 7 ?-band (Shapley Optical Survey; SOS; 
Mercuri o eral]|2006l ) and GALEX ultraviolet (|Rawle et all 
2008a). Inclusion in the high-S/N multi-object spectroscopy 
of iSmith. Lucev fc HudsonI (|2007D was a prerequisite for se- 
lection, and thirteen of the targets also have spectra from 
the N OAO Fundamental Plane Survey fNFPS: ISmith et all 
120041 ). 

The sample galaxies are all confirmed supercluster 
members, and were selected to probe a range of masses and 



age /metallicity combinations (as derived from previous cen- 
tral spectroscopy), and span three magnitudes in luminos- 
ity (m_R ~ 14 - 17 mag; Mr ~ -22.5 - -19.5 mag, using 
a distance modulus of 36.5). The galaxies all lie on the B- 
R colour-magnitude rel ation, and do not e xhibit unusually 
blue NUV-IR colours (|Rawle et al.ll2008al ). For consistent 
galaxy-to-galaxy comparisons, the PSF-corrected effective 
radius (ro) is a useful scale-length. We calculate from 
2MASS imaging, which is available for all our Shapley and 
Abell 3389 galaxies (see R08 for further details). The Shap- 
ley galaxies have effective radii of 1.7 - 9.7 arcsec (1.6 - 9.2 
kpc). Table [T] lists various parameters for the sample galax- 
ies. 

If present, nebular emission fills in the Balmer ab- 
sorption line s , cau sing H/3 to overestimate the galaxy age. 
ISmith et alj (|20o3) confirmed that Hq emission is a reli- 
able indicator of contaminating H/3 nebular emission. Cen- 
tral spectroscopy shows that eleven of our galaxies have no 
trace of Ha emission. Of the other eight, only three have 
EW(Ua) > 0.3 A (MGP012 9 MGP2083, MGP3 976), and 
none have EW{Ra) > 0.5 A. ISmith et al.l l|2009t) calculate 
that EW{}ia) = 0.5 corresponds to a maximum H/3 contam- 
ination of -0.1 A, which translates into an age offset of < 
15 per cent. At this low level, estimation of the H/3 emission 
would be unreliable, so we do not attempt corrections. 

T he targets were obser ved in April- July 2008 with VI- 
MOS l|LeFevre et al.1 12003| ) on the ESO VLT Melipal at 
Paranal. VIMOS was used in Integral Field Unit (IFU) mode 
with the high resolution, blue grism (dispersion = 0.51 A 
pixel"^), resulting in a wavelength range of 4150 - 6200 A, 
spectral resolution of 2.1 A (FWHM) and a field of view 
of 27 X 27 arcsec (0.67 x 0.67 arcsec per element). Obser- 
vations were made in dark time, with an average seeing of 
0.85 arcsec (FWHM). The hour-long observing blocks con- 
tained two 1022 s, spatially-offset exposures, and no sepa- 
rate sky frames. We aimed for a minimum S/N ~ 18 A-i 
at 5000 A for a 1 arcsec wide annulus at the effective radius 
(rc). To achieve this target for the five faintest objects (R 
< 16.0), the observing block was executed twice. A poten- 
tial issue with instrument fringing led to the repetition of 
several planned observations. Data quality was not in fact 
impacted, so all of the blocks have been included, increasing 
the co-added exposure time for some galaxies (detailed in 
Table©. 

2.2 Data reduction 

The initial data reduction was broadly similar to the tech- 
niques detailed in R08. However, there were some significant 
differences, the salient points of which are described below. 

Cosmic rays were removed from the raw frames 
using LACOSMIC (Ivan DokkumI l200ll ). Reduction to a 
sky-subtracted and wavelength calibrated datacube was 
achieved using the ESO common pipeline (EsoRex) and the 
IRAF routine imcombine. VIMOS is comprised of four de- 
tectors, so every exposure is split into four quadrants. Each 
quadrant was fiat-fielded and wavelength calibrated sepa- 
rately, using arc lamp exposures from immediately after the 
target observations. A master sky background spectrum was 
calculated for each quadrant by taking the median over all 
the spatial pixels in the outer region (rejecting bad pixels). 
This sky spectrum was then subtracted from each spatial 
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Table 1. Observed parameters of the Shapley Supercluster target galaxies. ID from lMetcalfe. Godwin &: PeachI l|l994l . MGP). Positions, 
central optical photometry (5 arcsec radius apertures) and morphology (by visual inspection) from SOS B- and i?-band imaging. J-band 
(total) magnitude and PSF-corrected circularised effective radius (ro) from 2MASS images. Ttot is the co-added VIMOS exposure time. 
Central S/N A"! ( averaged over the range 4700—5000 A) and central velocity dispersion (a km s ^; as measured in Section 12.31 1. both 
from an re/3 aperture. The final column indicates the number of radial bins used for the derivation of gradients (within r^), if the analysis 
was possible. 



ID 


RA 


Dec 




B-R 


mj 




Morph. 




S/N 






Bins 


(MGP) 


(J2000) 


(J2000) 


mag 


mag 


mag 


arcsec 




(s) 


A-i 


(km 


s 




# 








16 42 


1 '^7 


14 62 




Hi j DU 




14 1 


103.5 




U. 1 




MGP1189 


13:26:55.9 


-31:24:47 


15.46 


1.47 


13.50 


1.66 


E/SO 


4088 


44.4 


289.4 


± 


4.7 


4 


MGP1195 


13:26:56.0 


-31:25:27 


15.76 


1.42 


13.65 


3.25 


SO 


2044 


22.3 


115.5 


± 


3.1 


3 


MGP1211 


13:26:58.3 


-31:30:30 


16.83 


1.39 


15.19 


1.14 


E 


4088 


11.7 


88.9 


± 


9.6 




MGP1230 


13:27:00.3 


-31:21:05 


16.16 


1.30 


14.41 


2.07 


E/SO 


4088 


19.7 


57.5 


± 


3.8 


3 


MGP1440 


13:27:27.0 


-31:41:07 


15.92 


1.43 


13.97 


2.34 


E 


6132 


25.9 


187.6 


± 


5.5 


3 


MGP1490 


13:27:33.5 


-31:25:26 


16.40 


1.41 


14.24 


3.65 


E 


8176 


14.5 


76.8 


± 


4.8 


3 


MGP1600 


13:27:46.6 


-31:27:15 


15.02 


1.49 


12.96 


2.23 


SO 


2044 


43.5 


297.4 


± 


4.8 


5 


MGP1626 


13:27:48.5 


-31:28:46 


15.35 


1.55 


13.33 


1.66 


SO 


2044 


36.8 


259.5 


± 


6.8 


4 


MGP1835 


13:28:10.5 


-31:23:10 


15.19 


1.53 


12.93 


3.68 


E 


2044 


44.3 


202.1 


± 


4.3 


5 


MGP1988 


13:28:27.6 


-31:48:18 


15.76 


1.54 


13.69 


2.48 


E/SO 


2044 


26.9 


175.8 


± 


5.7 


4 


MGP2083 


13:28:38.7 


-31:20:49 


14.76 


1.65 


11.86 


9.62 


E 


4088 


47.8 


243.9 


± 


3.2 


7 


MGP2146 


13:28:47.2 


-31:41:35 


15.69 


1.42 


13.36 


4.61 


SO 


4088 


21.6 


121.3 


± 


4.0 




MGP2399 


13:29:14.5 


-31:39:11 


16.52 


1.39 


14.70 


2.33 


E 


6132 


11.8 


80.5 


± 


7.4 




MGP2437 


13:29:20.1 


-31:16:39 


15.74 


1.54 


13.84 


1.87 


E 


2044 


29.2 


186.7 


± 


6.1 


3 


MGP2440 


13:29:20.7 


-31:32:25 


14.99 


1.48 


12.38 


6.34 


E 


4088 


50.1 


227.0 


± 


3.5 


6 


MGP3971 


13:33:09.5 


-31:36:12 


15.77 


1.51 


13.84 


1.93 


E/SO 


2044 


21.0 


182.6 


± 


9.2 


3 


MGP3976 


13:33:11.7 


-31:40:09 


15.26 


1.52 


13.15 


2.52 


E/SO 


2044 


27.5 


252.6 


± 


5.3 


4 


MGP4358 


13:34:08.1 


-31:47:34 


16.86 


1.46 


15.11 


2.20 


E/SO 


6132 


9.0 


68.0 


± 


9.1 





element, pixel- by-pixel, after normalising to the total flux in 
the [01] A5577 A sky line. The final datacubes were produced 
by median combining the exposures in all of the observing 
blocks, followed by relative flux calibration with reference 
to a spectroscopic standard star. Approximate relative flux 
calibration is suflicient for a study of spectroscopic gradients 
in a sample of galaxies at a common redshift. 

2.3 Spatial binning and line index measurements 

As in iRawle et al.l (l2008lj ). the S/N per unit area is too 
low to attempt full two-dimensional mapping of the spectral 
features. Instead, a series of annuli are used to compress 
the data down to one spatial dimension (ellipticity-corrected 
radius) . 

The ellipticity and position angle of an ellipse fit to 
the broadband, wavelength-collapsed image of each galaxy 
formed the basis of the annuli used as radial bins. The annuli 
were arranged out to for each galaxy individually. Each 
annulus was required to have a S/N > 18 A~^ (average for 
the range 4700-5000 A) and a width greater than 0.7 arcsec. 
Central bins have a minimum semi-major axis of 0.5 arcsec. 
The galaxies were assigned between three and seven radial 
bins, depending on their effective radius and surface bright- 
ness. As an illustration, the adopted bins for four galaxies 
are shown in Figure [1] The average seeing disc size is simi- 
lar to the minimum annulus width, so generally each bin is 
largely independent of neighbouring bins. Spectra were also 
extracted from within three circular apertures, with radii 
of 1 arcsec (for direct comparison to previous central spec- 
troscopy; Section [H}, ro/3 (for physically consistent central 
parameters; Section [5} and (used in calculating general 
aperture corrections; Appendix 



The velocity dispersion, cr, was computed from the 
spectrum of each bin usin g the fitting routine ppxf 
(ICappeUari fc Emsellemll2004l ) and the lVazdekij (|l999l ) evo- 
lutionary stellar population templates. The templates cover 
a similar wavelength range to VIMOS (~4900 - 5500 A 
is used for fitting) and have a higher spectral resolution 
(FWHM = 1.8 A). Uncertainties on a were calculated using 
Monte Carlo simulations. 

We used INDEXI0 to measure atomic and molecular Lick 
indices from de-redshifted, flux-calibrated spectra at the VI- 
MOS instrument resolution (~50 km s~^). For measuring 
the gradients, model comparisons are performed in a rela- 
tive sense, so we do not calibrate to the Lick instrumental 
response. Uncertainties were calculated by indexf (as de- 
scribed in lCardiel et al.lll998l ) from error spectra estimated 
from the variance plane. We correct the measured line in- 
dices to account for both the instrument resolution and for 
broadening due to the observed velocity dispersion of each 
bin, adopting the same single-stage linear method used in 
R08. The technique avoids artificially broadening the ob- 
served spectra prior to index measurement, preserving any 
non-uniform noise properties. 



3 IFU CENTRAL SPECTROSCOPY 

One attraction of IFU data is the versatility in possible ap- 
proaches to the analysis. By extracting a spectrum from a 
circular aperture, it is trivial to mimic a single-aperture ob- 
servation. Using the 1 arcsec radius apertures introduced in 
Section 12.31 the VIMOS data can be compared directly to 

^ http:/ /www. ucm.es/info/Astrof/software/indexf/ 
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Figure 1. Adopted radial bins for four example galaxies: (left-to-right) MGP1600, MGP1835, MGP2440, MGP3976. Upper panels show 
broadband, wavelength-collapsed images from VIMOS, overlaid by the bin boundaries (green) and the effective radius (rc; red). Lower 
panels show the S/N (averaged over the region 4700 - 5000 A) within each annulus out to 10 arcsec; rc is indicated. 



two previous spectroscopic studies. For reference, these VI- 
MOS central apertures have an average S/N ~ 25 A"^ over 
the wavelength range 4700 - 5000 A. 

The first co mparison study is the NOAO Fundamen- 
tal Plane Survey l|Smith et al.ll2004 NFPS), which analysed 
1.5-2 hour observations from the Hydra multi-fibre spectro- 
graph on the 4 m CTIO telescope. The spectral resolution 
was ~3 A (FWHM) and the fibre diameter was 2 arcsec. 
Thirteen of the VIMOS target galaxies were included in 
NFPS, with an average S/N ~ 40 A'^ (5000 - 5500 A). 

The second is a deeper study from ISmith et al.1 (|2007h . 
which employed 8 hour observations from the AAfi fibre-fed 
multi-object spectrograph on the 3.9 m Anglo- Australian 
Telescope. The spectral resolution was ~3.5 A (FWHM) and 
the fibre diameter was again 2 arcsec. All of the VIMOS 
targets were observed by the AAfi study, with an average 
S/N ~ 100 A"^ (for 4700 - 5000 A). 



3.1 Central line strengths 

Figure [2] compares the VIMOS line strength measurements 
(for Fe4383, Fe4668, H/3 and Mgb5177) to those from the 
two studies introduced above, using the matched 1 arcsec 
radius aperture. Agreement between the VIMOS and NFPS 
values is generally very good, with a reduced ~ 1 (s-f" 
ter accounting for small mean offsets) . Measurements of the 
Fe4383 and H/3 lines agree well across all three instruments. 
Fe4668 and Mgb5177 exhibit comparatively large offsets be- 
tween AAQ. and both other studies, with the AAf2 results 
exhibiting weaker measured lines. 

For Fe4668, there is a (> 5(t) correlation between the 
VIMOS-AAn offset and the index value itself. The Fe4668 
index is one of the most susceptible to errors in the contin- 
uum level, as it is a particularly broad feature. Exploring fac- 
tors that could systematically affect the continuum level, we 
find no correlation between the offset and either airmass or 
position in the AAf2 multi-object spectrograph field. Aper- 
ture 'smearing' (the increase of the effective aperture size in 



0.11 21.67 / 18 



0.22 8.30 / 12 





3 4 5 6 
Fe4383 (NFPS) 



3 4 5 5 
Fe4383 (Mil) 



O.SO 34.83 / 18 



3 4 5 6 
Fe4383 (AAtl) 



1.10 19.68 / 12 















X . , . 





4 6 8 
Fe4668 (NFPS) 



4 6 8 

Fe4668 (AAQ) 



0.01 21.95 / 



Fe4568 (AAQ) 
-0.02 11.03 / 12 




1 1.5 2 2.5 
H/J (NFPS) 


1 1.5 2 2.5 
H/! (AAtl) 


0,06 


15.20 / 12 


0,30 27.14 / 18 










5- ■ 


y. : 



1 1.5 2 2.5 
H/S (AAB) 



0.25 20.74 / 12 




Mgb5177 (NFPS) 



Mgb5177 (AAD) 



Mgb5177 (AASl) 



Figure 2. Comparison of VIMOS central spectroscopy to NFPS 
(left column) and AASl (central column). Plots of NFPS versus 
AAf! index measurements for these objects are displayed in the 
right column. Dashed lines are equality; solid lines shows the offset 
with the minimum x^- The values above each panel show the 
offset and minimised / degrees of freedom. 
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offset = -0.03 offset = +0.04 offset = -0.03 




central log(age/Gyr) (basic) central [Z/H] (basic) central [a/Fe] (basic) 

Figure 3. Comparison between VIMOS spectroscopic central age (left panel), metallicity (centre panel) and a-abundance (right panel) 
derived using two different index sets: basic (H7F + H/3 + Fe4383 + Fe5015 + Mgb5177; top line in Table EJ; master (basic + USF + 
Fe4531 + Fe4668 + Fe5406; bold line in Table [2j. Solid line shows the weighted mean offset; dashed line is equality. The mean offset is 
displayed above each panel. 



Table 2. Mean values of the central age, metallicity [Z/If] and o-element abundance [o/Fe] derived using different index sets in the 
model inversion. The mean central stellar population parameters derived from AAf2 observations for the same 19 galaxies, are displayed 
in the final row for comparison. 



Index set 


(log(age/Gyr) ) 


( [Z/H] > 


( ["/Fe] > 


H7F 


f H/3 + Fe4383 + Fe5015 + Mgb5177 (basic) 


0.97 ± 0.04 


0.12 ± 0.04 


0.24 ± 0.03 


basic 


+ Fe4531 


0.97 ± 0.04 


0.12 ± 0.04 


0.24 ± 0.03 


basic 


+ Fe5406 


0.94 ± 0.04 


0.14 ± 0.04 


0.23 ± 0.03 


basic 


+ Fe4668 


0.95 ± 0.04 


0.14 ± 0.04 


0.24 ± 0.03 


basic 


+ H<5F 


0.97 ± 0.04 


0.11 ± 0.04 


0.24 ± 0.03 


basic 


+ Fe4531 + Fe5406 


0.95 ± 0.04 


0.14 ± 0.04 


0.22 ± 0.03 


basic 


+ Fe4531 + Fe4668 


0.95 ± 0.04 


0.15 ± 0.04 


0.24 ± 0.03 


basic 


+ Fe4668 + Fe5406 


0.94 ± 0.04 


0.15 ± 0.04 


0.23 ± 0.03 


basic 


+ H(5F + Fe4531 + Fe5406 


0.93 ± 0.04 


0.15 ± 0.04 


0.22 ± 0.03 


basic 


+ H(5F + Fe4531 + Fe4668 


0.94 ± 0.04 


0.14 ± 0.04 


0.23 ± 0.03 


basic 


+ H(5F + Fe4668 + Fe5406 


0.93 ± 0.04 


0.16 ± 0.04 


0.22 ± 0.03 


basic 


+ Fe4531 + Fe4668 + Fe5406 


0.94 ± 0.04 


0.16 ± 0.04 


0.22 ± 0.03 


basic 


+ H<5F + Fe4531 + Fe4668 + Fe5406 (master) 


0.93 ± 0.04 


0.16 ± 0.04 


0.22 ± 0.03 


AAn [VIMOS targets] basic + HSF 


0.96 ± 0.04 


0.10 ± 0.03 


0.24 ± 0.03 



co-added multi-object observations due to fibre placement 
or telescope pointing errors) may explain a portion of the 
offset, but cannot account for the entire correlation. While 
the responsible factor has not been identified, the agreement 
between VIMOS and NFPS measurements, and their similar 
offsets from AAfi, suggests that the problem may lie with 
the multi-object spectra from AAfi. 

3.2 Central stellar populations 

Stellar populations, characterised by age, total metallicity 
(Z/H) and a-element over-abundance (a/Fe), were derived 
through comparison of the measured line indices to predic- 
tions from single burst evolution models. We use a grid inver- 
sion t echni que similar to that of IProctor, Forbes fc BeaslevI 
l|2004l ) and lSmith. Lucev fc HudsonI l|2009l ). This involves a 
non-linear minimisation over a set of indices, using a 
cubic-spline interpolation in age-[Z/H]-[a/Fe], and a lin- 
ear extrapolation for points l ying outside the model grids. 
We co mpare to the models oflThomas. Maraston fc Bended 
l|2003l ) and iThomas. Maraston fc Koriil (|2004l ). hereafter 
TMBK. 



Nine indices were used in the inversion (the 'master' 
set): three Balmer lines (H(5F, H7F, H/3), five iron lines 
(Fe4383, Fe4531, Fe4668, Fe5015, Fe5406) and the a/Fe sen- 
sitive Mgb5177. H(5A and H7A are considered redundant as 
they mostly duplicate the response of H5F and H7F. Iron 
lines towards the red end of the VIMOS wavelength range 
(Fe5270, Fe5335, Fe5709, Fe5782) are excluded due to con- 
tamination by sky features, such as the [OI] A5577 line. 

To test the stability of the inversion technique, we de- 
rived central stellar population parameters using several dif- 
ferent sets of indices. Starting from a 'basic' set of five indices 
(H7F, H/3, Fe4383, Fe5015, Mgb), thirteen sets were created 
by adding combinations of other lines in the 'master' list. 
The mean characteristics of the 19 galaxies are robust with 
respect to the different index sets (Table [2} . The mean age 
of the sample varies by less than 5 per cent, and the mean 
metallicity and a-abundance deviate by only 0.04 dex and 
0.02 dex respectively. The mean central parameters of the 
same galaxies derived from AAfl spectroscopy (final row. 
Tabled are well matched to the VIMOS values. Individu- 
ally, the derived central stellar populations of galaxies are 
also independent of the index set used. Figure [3] compares 
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Figure 4. Comparison between derived stellar population parameters from VIMOS (using master index set) and AAQ. Layout as in 
Figure |3] The distributions and mean values in each parameter are well matched between the two observations, although the errors 
appear slightly underestimated, particularly in [Z/H]. 



values computed from grid inversions using the basic and 
master sets. 

For the remainder of this study, we use the master set 
of indices (bold row in Table [2} to calculate stellar popu- 
lation parameters. Figure |3] compares the central values of 
age, metallicity and a-abundance derived from VIMOS to 
the equivalent values from AAfi. Although the mean value 
for each parameter is well matched between observations, 
the reduced x'^ is ~ 4 for age and [a/Fe], and ~ 6 for metal- 
licity. It appears that the formal errors, propagated through 
the index grid inversion by calculating the full error ellip- 
soids from the surface, are underestimated. Errors on the 
indices themselves are well determined, as noted in Section 
13.11 (Figure [2} . Consequently, in the following sections we 
quote the uncertainty on stellar population parameters as 
derived from Monte Carlo simulations of the index values. 



4 GRADIENTS 

4.1 Line strength gradients 

Line strength gradients were measured using an unweighted 
linear least-squares fit to all the elliptical bins within the 
effective radius (re, see Section [2.3p . such that VI = AI / 
Alog(r/rc), where I is any line strength index. Most previous 
studies of radial trends have used the log(r/ro) convention, 
resulting in straight forward comparisons. A small fraction of 
galaxies would be better fit in linear radial space; a weighted 
fit to these in log(r/re) would bias the overall gradient to- 
wards the radial trend of the inner region, where the index 
values are better determined. The unweighted fit ensures a 
meaningful parametrization of the average gradient out to 
Tc in all galaxies. 

Five galaxies (MGP0129, MGP1211, MGP2146, 
MGP2399, MGP4258) have S/N < 18 A"^ at r^, due to 
their lower surface brightness. These galaxies are excluded 
from further analysis as the quality of the spectra is 
insufficient to calculate reliable gradients. 

Generally, the metal lines show negative gradients, i.e. 
weaker lines at larger radii, while the Balmer lines have small 
positive gradients (Figure [S}. The gradients do not exhibit 
simple trends with the central velocity dispersion, although 
the exclusion of most of the small galaxies, due to low surface 



brightness, results in a poor sampling of the low mass end 
(log a < 2.2). 



4.2 Stellar population gradients 

For each of the 14 galaxies with measured gradients, the lin- 
ear fits described above were extrapolated to estimate the 
index values at r = O.lxrc and r = Vc- The correspond- 
ing stellar population parameters for these two points were 
derived using the grid inversion technique described in Sec- 
tion 13.21 and the 'master' index set. Finally, gradients in 
age, metallicity [Z/H] and a-element enhancement [a/Fe] 
were calculated from the difference between these inner and 
outer points. 

This process differs from the approach used in R08, 
where stellar population parameters were calculated for each 
annulus, and gradients were fit in a similar manner to the 
line indices. To test the sensitivity of the results to the choice 
of technique, Monte Carlo simulations were computed for 
the measured line strengths (based on the indexf errors), 
and stellar population gradients were derived for each galaxy 
using both methods. The new method proved to be more 
stable, relying on 2-3 times fewer model inversions and not 
requiring a linear fit to noisy stellar population values. The 
average simulated gradient for each parameter is method- 
independent. Errors quoted for the stellar population gradi- 
ents are computed from these Monte Carlo simulations. 

The point spread function (PSF) of the observations 
may lead to an underestimation of the gradients in the small- 
est galaxies, as central light is smoothed into neighbour- 
ing bins. To quantify this effect, we constructed a simple 
two dimensional model with a known metallicity gradient 
convolved with an r'^^* luminosity profile. Seeing conditions 
were simulated by smoothing the model with Gaussian fil- 
ters of widths in the range of the measured PSFs (mean 
FWHM = 0.85 arcsec; maximum FWHM = 1.2 arcsec). For 
each smoothed model, the 'observed' gradients were derived 
from radial bins similar to those used for the real data. In 
the worst realistic case, with PSF FWHM — O.Qxro, the 
measured gradient was fiatter than the actual gradient by 
25 per cent. For the majority of the Shapley observations, 
the effect is estimated to be at the 5-10 per cent level, so 
no corrections have been applied. 
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Figure 5. Absorption line strengtli gradients versus central logcr (rc/3 radius aperture). The solid line represents the median gradient 
(also displayed at the top of each panel); dashed line indicates zero gradient. 



On average, the galaxies in the Shapley sample have 
a negative gradient in all three stellar population parame- 
ters. The mean age gradient is -0.09 ± 0.07, although there 
are individual galaxies with significant gradients in either 
direction. Both metallicity and [a/Fe] have marginally sig- 
nificant mean negative gradients; { V[Z/H] ) = -0.11 ± 0.05 
and {V[Q/Fe]) = -0.06 ± 0.03. The top row of Figure [g] 
presents the distributions of the gradients for the 14 Shap- 
ley galaxies and Table lB2l in Appendix |B] lists their values. 



4.3 Trends in gradients 

This section presents the correlations between the stellar 
population gradients and central parameters. The latter are 
measured from the ro/3 apertures, as listed in Table IBll 
Physical interpretation of the observed dependencies is de- 
ferred to Section [5l 

The Shapley results will be comb ined with gradients 
derived for the galaxies in Abell 3389 (|Rawle et al.ll2008bl . 
R08). The raw observations from R08 have been re-analysed 
using the current methodology, producing stellar population 
gradients for 11 of the 12 galaxies (D53 has insufficient S/N 
at Te). None of the general conclusions in R08 are altered by 
the re- worked analysis. The second row of Figure |6] presents 
the re-calculated gradient distributions for the R08 galaxy 
sample (hereafter, A3389), and the stellar population char- 
acteristics are also presented in Appendix [C] The combined 
set of galaxies from Shapley and A3389 (termed VIMOS) 
have {Vlog(age)) = -0.02 ± 0.06, {V[Z/H]) = -0.13 ± 
0.04 and ( V[a/Fe] ) = -0.10 ± 0.04. 

We explore the existence of relationships between pairs 
of central parameters and gradients, using a simple linear 



regression. The gradient is treated as the dependent vari- 
able, as initially our interest lies in determining how well 
the gradients can be predicted from knowledge of the cen- 
tral properties. Fits were calculated without weighting for 
errors as intrinsic scatter may be significant. Most pairs of 
parameters do not show a significant trend in the VIMOS 
data. TableOlists the slope and significance of each. For ex- 
ample, none of the gradients have a simple dependency on 
central velocity dispersion. (We explore the V[Z/H] versus 
central loga relation further in Section 15.21 see Figure [T]) 
Two significant correlations stand out; the metallicity gradi- 
ent is related to the central metallicity (central upper panel, 
Figure[S| with a slope of -0.63 ± 0.20 (significance of 3.1 a); 
the age gradient is also correlated with the central metal- 
licity (upper panel. Figure |9]), with a slope of 1.09 ± 0.28 
(3.9 o"). Conversely, neither the age gradient nor the metal- 
licity gradient depend significantly on central age (upper-left 
panel Figure |8] and upper panel Figure [TOl respectively) . 

The final row of Table |3] refers to the relation between 
age gradient and metallicity gradient. In this case, it is un- 
clear which parameter is dependent on the other. Hence, 
the trend is calculated from the bisecting linear relation of 
the two simple regressions (age on metallicity; metallicity on 
age) . A strong relation is found in the VIMOS sample, with 
a slope of -0.69 ± 0.10 (significance of 7.0 a; upper panel. 
Figure [TT]). 



5 DISCUSSION 

Many factors potentially contribute to the observed trends 
in stellar population gradients. Hence, there is a risk of 
reaching a misleading conclusion due to sparse sampling of 
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Table 3. Summary of the relations between central parameters (within re/3 apertures) and gradients shown in Figures [Tl- llOl and/or 
mentioned in the text. The slope, error on the slope and significance (in parenthesis) are calculated by minimising the residuals in the 
Y-direction. The columns display trends for A3389, Shaplcy, VIMOS (= A3389 + Shapley) and the comparison studies (M03 and SB07). 
Significant trends in the VIMOS sample (> 3o-) are highlighted in bold type. The final row shows the slope for the bisector of the two 
simple regressions between age and metallicity gradients (Figure [TTJ. 
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Figure 6. Histograms showing the fractional distribution of 
gradients in (loft-to-right) log(age/Gyr), [Z/H] and [o/Fe] for 
galaxies in each sample: (top-to-bottom) Shapley, A3389 (re- 
analysed R08 data), M03, SB07 (see Section |ETJ. Mean gradi- 
ents for each sample are marked by dashed lines, and are shown 
in the top-left of each panel. The mean gradients for the VI- 
MOS (Shapley+A3389) sample are marked by a dotted line in 
the top two rows, and presented at the top of each column. Back- 
ground shaded histograms are for the overall distribution (Shap- 
ley+A3389+M03+SB07) . 



a complex parameter space. To minimise this possibility, we 
test the validity of the trends found in Section [4.31 by com- 
paring to data from two published early-type galaxy studies 
(detailed in Section [571} . In the first instance, the relation- 
ship between central velocity dispersion and metallicity gra- 



dient is explored (Section l5.2p . We then assess the predictive 
power of the central stellar population parameters in esti- 
mating the gradients (Section 15. 3p . Section [5.41 investigates 
the interplay between the gradients in age and metallicity. 

5.1 Comparison studies 



The 
are 



two main c ompa rison studies considered 
Mehlert et al.l (|2003| . hereafter M03) and 



ISanchez-Blazguez et"al] (|2007| . SB07). Both used long 
slit spectroscopy to derive stellar population gradients, 
emp loying similar t echniq ues and models to our analysis. 

iMehlert et al.l (120031 ) derived stellar population gradi- 
ents for 32 early-type galaxies in the Coma cluster. The 
authors rebinned the spectra to achieve a S/N > 30 
per bin out to r^, and computed the stellar population 
parameters using an index-pair iteration technique, with 
H/3, Mgb5177 and (Fe) = (Fe5270+Fe5335)/2, in conjunc- 
tion with the TMBK models. Their original analysis ex- 
cluded the two cD galaxies (NGC4889, NGC4874) and an 
SO (NGC4865), as they lie beyond the model grids, and 
NGC4876 due to systematically larger errors. We follow suit, 
leaving 28 galaxies in the sample. M03 defined central val- 
ues as the average along the major axis inside O.lro/Vl — 
where e is the ellipticity. The distributions for the gradients 
are shown in the third row o f Figu re [6] 

ISanchez-Blazguez et al.l (|2007l l studied a sample of 11 
local early-type galaxies in various environments (field, 
group, Virgo cluster), rebinning the spectra to achieve a S/N 
> 50 A^^ out to 2rc. The authors derived stellar population 
gradients for each galaxy using a minimisation technique, 
with 19 Lick indices and the TMBK models. NGC2865, a 
merger remnant, contains an extremely young core, and has 
been excluded from our analysis. For central values, SB07 
used an aperture size of 1.5 arcsec x ro/8. The final row 
of Figure [6] displays the distribution of gradients for the 10 
SB07 galaxies. 

In general, the average gradients of the four studies are 
very similar. The Shapley sample has a more steeply nega- 
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tive mean age gradient, and the combined VIMOS sample 
has a larger mean V[a/Fe] than either of the long-slit in- 
vestigations. The mean negative metallicity gradient from 
SB07 is nearly twice as steep as in the other samples. The 
authors derive large negative V[Z/H] gradients for galaxies 
in both high and low density environments, so the difference 
is not due to their inclusion of non-cluster targets. 

5.2 Metallicity gradient versus central log a 

The combined VIMOS sample does not show a significant 
correlation between metallicity gradient and the central ve- 
locity dispersion (upper panel. Figure [7|. The lack of a sim- 
ple linear relationship is confirmed in the comparison stud- 
ies. As we discuss below, evidence has been previously re- 
ported in favour of the [Z/H] gradient depending on mass 
(via various observable proxies). 

A simple model of galaxy formation asserts that en- 
richment gradients depend on galactic wind efficiency, 
which scales wi t h th e potential and hence galaxy 
mass (iK obavaslu 2004 ) . I n a study of 42 early- types, 
ICarolloTPanziger fc BusonI (|l993l ) measured the gradient in 
the Mg2 absorption line (a simple proxy for metallicity) and 
discovered that it became shallower with decreasing cen- 
tral velocity dispersion (galaxy mass) for low-mass systems 
(log cr < 2.2 6). High-mass galaxies, howev er, showed no such 
correlation. iKobavashi fc Arimotol (|l999l ) stu died the Mg2 
gradients of 80 galaxies (above and below the lCaroUo et ahl 
transition point), but were unable to find a statistically sig- 
nificant trend with any mass estimator (central velocity dis- 
persion, B-band luminosity, effective radius or the dynami- 
cal mass). 

In contrast, iForbes. Sanchez-Blazguez fc Proctoij 
l|2005ll claimed a significant relation between metallicity 
gradients and galaxy mass (if-band luminosity, central 
velocity dispersion and dynamical mass) fo r a literature 
sampl e of Coma cluster galaxies. Recently, ISpolaor et al.l 
(|2009l ) have resurrected the idea of a mass transition point 
after observing 14 low mass early-types in nearby clusters 
(Virgo and Fornax). These authors found that for low mass 
galaxies (log cr < 2.15), metallicity gradients are remarkably 
tightly correlated with velocity dispersion, in the sense that 
(negative) V[Z/H] becomes shallower with decreasing logcr 
and positive at very low logcr. They conclude that while 
some higher mass galaxies have steep negative V[Z/H], 
there is also a visible downturn and a markedly broader 
scatter. Physically, the transition point may mark the mass 
above which galaxy formation histories are dominated by 
late-epoch , gas-poor merging , which dilute the metallicity 
gradient (jHopkins et all [20090 . The magnitude of the 
dilution is related to the individual merge r properties and 
the m etallicity gradients of the progenitors l|di Matteo et al.l 
I2OO9I ). causing the observed inflated scatter at larger masses. 

As already noted. Figure [7] confirms the large scatter 
and apparent absence of correlation for high mass galaxies. 
Unfortunately, the A33 89, M03 and S B07 studies do not in- 
clude galaxies below the lSpolaor et afl mass transition point. 
Low mass galaxies are also under-represented in the Shap- 
ley sample, as several small galaxies have been excluded by 
low surface brightness at re. However, the remaining galax- 
ies with logcr < 2.15 are not inconsi stent with the mass - 
metallicity gradient trend reported in ISpolaor et all l|2009f) . 
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Figure 7. Metallicity gradient versus central velocity dispersion 
(logcr). Upper panel displays the VIMOS data (Shapley = black 
circles; A3389 = red squares). Lower panel shows VIMOS data 
together w ith M03 (blue open triangles), SB07 (magenta open 
circles) and lSpolaor et al] l|2009, S09; green open stars - measured 
from their figure 1). The two low mass galaxies from Shapley are 
consistent with the tight relation for low mass early-types in S09. 



5.3 Gradient trends with central values 

We now explore the relation between central stellar popu- 
lation parameters and gradients. The VIMOS sample shows 
no significant correlation of the age gradient with central 
age (upper left panel of Figure [8]), or in the analogous rela- 
tion for Q!-enhancement (upper right panel of Figure |8]). The 
comparison studies confirm the lack of trends in these cases 
(lower panels), as shown in the final two columns of Table 

m 

In contrast, the metallicity gradient in VIMOS galax- 
ies is related to the central [Z/H], with a slope of -0.63 ± 
0.20 (3.1 cr significance; upper central panel, Figure|8]|. This 
relation is confirmed in the M03 sample, but is not signif- 
icant in SB07. A sim ilar relation, with a slope — -0.71 ± 
0.34, was reported bv lMorelli et al.1 l|2008l ) for the bulges of 
nearby spiral galaxies, which are s ometimes considered aki n 
to low-luminosity ellipticals (e.g. iThomas" fc Daviesll2006h . 
The trend intercepts the origin, and the orientation of the 
gradients suggest that in general, galaxies tend towards so- 
lar metallicities ([Z/H] = 0) at ~ 2-3 r^. In rapid dissi- 
pative collapse, in-falling gas becomes enriched by evolved 
stars, and contributes to the metal-rich star s forming in the 
core. (Recent models by iPipino et al. I I2OO8I predict that as 
much as 90 per cent of metals in central stars originated in 
outer regions.) The onset of galactic wind depends on lo- 
cal escape velocity and occurs later in inner regions. As a 
result, star formation abates earlie r in outer regions caus- 
ing lower [Z /H] at re than the core l|Arimoto fc Yoshiilll987l : 
lMatteuccilll994l L Both of these processes work to strengthen 
the metallicity gradient and raise the central Z/H. 

The age gradient is correlated with central metallicity 
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Figure 8. Stellar population gradients versus the corresponding central value for (left-to-right) age, metallicity and o-abundance. Upper 
panels present the VIMOS data (Shapley = black circles; A3389 = red squares), together with best fit lines (dashed; solid black = 
combined VIMOS) as described in Section 14.31 Lower panels show results from M03 (blue triangles) and SB07 (magenta circles). Dash- 
dot line indicates the best fit to the combined comparison sample (M03 -I- SB07); solid line displays the combined VIMOS best fit. The 
slope and significance of the best fits are displayed to the top-loft of the panels, and in Table [3] 



in the VIMOS sample, with a slope of 1.09 ± 0.28 (3.9 a 
significance). The relation is such that galaxies with more 
metal rich centres have stronger age gradients. However, the 
relation is not significant in either comparison study, or in 
the A3389 sample alone. The metallicity gradient versus the 
central age also displays a marginal trend, in which centrally 
young galaxies favour stronger negative metallicity gradi- 
ents, and galaxies with old cores exhibit weaker metallicity 
gradients. 

To summarise, the metallicity gradient has a strong de- 
pendence on the central metallicity. Although further corre- 
lations between parameters are only marginal, some tenta- 
tive generalisations can be noted. On average, galaxies with 
metal rich cores have steeper negative metallicity gradients 
and Vlog(age) > 0. The observed trends are such that galax- 
ies with solar (or sub-solar) metallicity centres tend towards 
flat metallicity and negative age gradients. These results 
suggest, or are possibly driven by, a mutual correlation be- 
tween the age and metallicity gradients. 

5.4 Gradient— gradient relation 

The gradients in age and metallicity are strongly anti- 
correlated in the VIMOS sample, with a slope = -0.69 ± 
0.10 (we discuss the issue of correlated measurement er- 
rors in the next paragraph). The relation is confirmed by 
both comparison studies (Figure 1111 and the final row of 
Table |3]) . The trend is reminiscent of the central parame- 
ter age-met allicity relation, which also has a strong anti- 
correlati on with a slope of ~0.7 and a significant intrinsic 
scatter (|j0rgensenlll999l ). However, the central relation is 
a two-dimensional projection of t he three paramete r age- 
mass-met allicity plane CZ-plane': iTrager et ahlboOOl ). with 
much of the scatter originating from the mass dependence. 
In contrast, the trend in the gradients has no observed de- 
pendence on velocity dispersion. 

The strong ant i- correlation between Vlog(age) and 




U03 : 
SB07 : 
M03+SB07 : 



22 ± 0.21 (I.Otj) 
20 ± 0.12 h.7G) 
24 ± 0.13 (1-9(^) 



AM03 

OSB07 




0.2 0.4 
central [Z/H] 

Figure 9. Age gradient versus central metallicity. Layout as in 
each column of Figure [8] 



V[Z/H] suggests that generally, a positive age gradient is 
coupled with a strong negative metallicity gradient, and a 
negative age gradient is usually associated with a weak or 
flat metallicity gradient. However, the use of multiple ab- 
sorption line indices to derive the age and metallicity does 
not entirely erase the de generacy, with a co rrelation domi- 
nating the errors (see iKuntschner et al.ll200ll for further dis- 
cussion) . It is prudent to calculate whether the strong trend 
seen in Figure [11] could be caused by correlated measure- 
ment errors. 

The upper panel of Figure [12] again presents the gradi- 
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0.4 0.6 0.8 1 1.2 1.4 

central log(age/Gyr) 

Figure 10. Metallicity gradient versus central age. Layout as in 
each column of Figure |H] 




Vlog(age/Gyr) 

Figure 11. Relation between the age and metallicity gradients. 
Layout as in each column of Figure \E\ 



ents for the Shapley and A3389 samples. On this occasion, 
the simple one-dimension (Ict) error bars in each axis have 
been replaced by the 1 a error ellipses, based on the distri- 
bution of the Monte Carlo realizations. The error ellipses 
are clearly orientated in the direction of the correlation. 

The lower panels of Figure [T2] show the distribution of 
the galaxies projected onto the best-fitting line from the up- 
per panel. We consider a null hypothesis in which the scatter 
around the zero value is entirely due to the correlated errors. 
The Shapley sample does not disprove this hypothesis, with 
an estimated probability of g = 0.2. For the A3389 sample in 




0.5 

Vlog(age/Gyr) 




-0.5 0.5 

Projected onto the line of best fit 

Figure 12. Upper panel: The age gradient-metallicity gradient 
relation for Shapley (solid points) and A3389 (open points), with 
the error estimates represented by ellipses (Shapley = dark grey; 
A3389 = light grey). The best fit line is shown. Lower panels: For 
each sample (top-to-bottom: Shapley, A3389, VIMOS = Shap- 
ley -I- A3389) the points are projected onto the best fit line and 
shown as a displacement from the median value (y-axis position 
is arbitrary). The one-dimensional chi-squared, degrees of free- 
dom and q, the probability that the distribution could be entirely 
attributed to the effect of correlated errors, are displayed. 



iRawle et al.l l|2008bl ). we concluded that the null hypothesis 
was incompatible with the data, i.e. that although there is 
a sizeable correlation of errors, there is also a real intrinsic 
scatter in the same direction. The newly derived age and 
metallicity gradients for the A3389 sample confirm this con- 
clusion, as does the combined VIMOS sample. While not 
contrary to the null hypothesis, the Shapley sample is also 
not inconsistent with the presence of intrinsic scatter. To 
explain the age-metallicity gradient relation for the VIMOS 
sample without invoking intrinsic scatter, the errors in the 
stellar population gradients would have to be a factor of ~3 
larger. Such an underestimation is unlikely, as the gradient 
errors, derived via Monte Carlo simulations, scale approxi- 
mately linearly with the well determined index uncertainties. 
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1.3 1.4 1.5 1.6 1.7 1.8 
from imaging 

Figure 13. Comparison of B-R colours from imaging with those 
derived from the VIMOS stellar population parameters. Each 
point corresponds to a matched aperture from a target: ro/3 ra- 
dius aperture (black), radius aperture (red), innermost annulus 
from VIMOS radial binning scheme (green), other annuli from VI- 
MOS radial binning scheme (blue). Equality is shown by the solid 
line. The red {B-R > 1.75) predictions from VIMOS are all from 
MGP2437. 



5.5 Colour gradients 

Many studies of broadband colour gr adients have ignored 
the possibility of a ge gradients (e.g. IPeletier et al.l [l990l : 
iTamura et al.l 1 20001 ). but here we find that Vlog(age) is 
often inconsistent with zero. By measuring colour gradi- 
ents in multiple bands, recent studies have attempted to 
incl ude the age grad ient. From SDSS and 2MASS photome- 
try, ImiEalJI^oi) derived average V[Z/H] = -0.25 ± 0.03 
and Vlog(age) = 0.02 ± 0.04. However, the age and metal- 
licity predictions are strongly correlated due to the usual 
degeneracy in broadband colours. For instance, a negative 
B-R colour gradient could be the result of a negative gra- 
dient in either metallicity or age. Additionally, our study 
has found that age and metallicity gradients are intrinsically 
anti-correlated, affecting the colour gradient in opposite di- 
rections. 

From a known stellar population and given models, the 
broadband colour can be predicted. In this section, we com- 
pare B-R colour gradient predictions, calculated from VI- 
MOS spectroscopic age and metallicity, to those observed 
in B- and f?-band im ages from the Shapley Optical Survey 
l|Mercurio et al.l |2006| . SOS). The pixel size is 0.24 arcsec, 
and the average seeing was ~0.8 arcsec (FWHM) in the 
_B-band and ~0.7 arcsec for R, corresponding well to the 
spatial information from VIMOS (0.68 arcsec sampling and 
an average seeing of 0.85 arcsec). For consistency, the B-R 
gradients were computed using a method as similar as pos- 
sible to the calculation of VIMOS line strength gradients. 
B- and i?-band magnitudes were measured for the same el- 
liptical annuli and circular central apertures (r-c/3 and r^) 
described in Section 12.31 Raw photometry was calibrated 




-0.2 -0.1 0.1 
from imaging 

Figure 14. Comparison between the B-R colour gradients from 
imaging and VIMOS stellar population parameters. Equality is 
indicated by the solid line. The arrow disp l ays th e mean \7B-R. 
= -0.08 ± 0.01 mag dex"! from lWu et"ai1 ll2005ll . 

onto the Johnson (Vega) system using the SOS image ze- 
ropoints, and c orrected for galactic extinction with the red- 
dening maps of ISchleeel. Finkbeiner fc David ()l998[ fl. 

B-R colours were predicted from the VIMOS spec- 
troscopic age and metall icity, using the e volutionary stel- 
lar population models of iMarastonI l|2005l . the basis of the 
TMBK models), which assumes the Kroupa IMF and so- 
lar Q-abundance. These predictions (for individual apertures 
and annuli) agree well with the photometry, with an overall 
rms dispersion of 0.068 mag (Figure [13)) . The innermost bins 
have been distinguished from other annuli, as they are most 
likely to be affected by the PSF (aperture diameter « PSF 
FWHM)Q 

The best VIMOS predictions are for the apertures 
(red points in Figure [T3}, which have high S/N and are not 
affected by seeing. The comparison has an rms dispersion 
— 0.056 mag, which is decreased to 0.039 mag if the sin- 
gle outlier (MGP2083: prediction = observed - 0.16 mag) 
is excluded. This scatter is consistent with that found by 
ISmith et al.l (|2009l ). who studied ~230 galaxies in Shapley 
using AASl spectroscopy and SOS aperture colours. To ac- 
count for the effect of a-abundance on the colour, that study 
empirically derived the relation SB-R — -0.14 x (5[Q/Fe]. 
Using this, our scatter for the apertures is further reduced 
to 0.034 mag. 

The three outlying points with very red predicted 
colours from VIMOS {B-R > 1.75) are all from MGP2437 
(rc/3 aperture and the two innermost annuli). The red 
colour prediction results from a large spectroscopic cen- 
tral age estimation (log(age/Gyr) = 1.412; lying beyond the 

^ As=4.315xE(B-y), Aji=2.e73x E{B-V) 

* The innermost 'annuli' are the central ellipses adopted during 
the spatial binning described in Section l2.3l In several cases they 
are very similar to the circular rc/3 central apertures. 
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model grids). The extreme age is not observed in the larger 
1 arcsec radius aperture, so the cause must be localised in 
the core. The anomaly appears to be due to extremely weak 
H(SF and H7F measurements, most likely resulting from a 
combination of poor S/N at the blue end of the spectrum 
and an unusually deep G4300 line (which lies in the blue 
'continuum' sideband of H7F). 

Colour gradients (for both observed and predicted 
colours) were calculated using an unweighted linear fit to 
all bins within the efi'ective radius: VB-R = AB-R / 
Alog(r/ro). We compared the predicted colour gradients 
to those observed, as shown in Figure 1141 The predictions 
are reasonably consistent with the observations, with a re- 
duced ~ 3.4 (rms scatter of 0.08 mag). Accou nting for the 
Vfa/Fe] efi'ect on the colour gradient (using the lSmith et al.l 
formula above) has only a marginal efi'ect, as the mean cor- 
rection is < 0.01 mag. The sample mean is well matched to 
themean VB-R = -0.08 ± 0.01 mag dex~^ from lWu etHI 
(|2005f ). 



5.6 Towards a physical explanation 

Drawing together the observed trends, we propose the fol- 
lowing tentative explanation. Galaxy mass dictates both 
the ef ficiency of gas in- fall during rapid dissipative col- 
lapse l|Pipino et ahllioosl ). and the delayed onse t of galac- 
tic w i nds in the core c ompared to outer regions l|Matteuccil 
1 19941 : lKobavashill2004h . The former process raises the cen- 
tral metallicity, while the latter inhibits metal enrichment at 
larger radii. Together, they ensure that the strength of the 
negative metallicity gradient correlates with central metal- 
licity, and that both parameters generally increase with ve- 
locity dispersion. For low mass galaxies {a < 140 km s~^) 
these are the pre vailing processes, dr i ving the relation fo r 
loga < 2.15 (e.g. ICaroUo et al.lll993l : ISpolaor et al.ll2009^ . 
Above this transition point, mergers play a dominant role. 

Gas-rich merging may prompt a central starburst, re- 
sulting in an increased ce ntral metallicity a nd a more nega- 
tive metallicity gradient l|Kobavashill2004h . reinforcing the 
strong V[Z/H] -central [Z/H] relation. The starburst also 
creates a positive age gradient (M ihos fc Ifernquist 19 94). 
which may be the origin of the Vlog(age)-central [Z/H] re- 
lation. However, the logarithmic age gradient from a cen- 
tral starburst is diluted over time, introducing a depen- 
dency on the epoch of the last merger, and weakening the 
Vlog(age) correlations. Late-ep och, gas-poor mergers tend 
to fla tten metallicity gradients ljWhitelll980l : iHopkins et al.l 
I2OO9I ). with the exten t of the dilution d ependent on the prop- 
erties of the merger l|di Matteo et al. 2009). This increased 
scatter in V[Z/H] weakens all of the correlations further. 



6 CONCLUSIONS 

The internal gradients of early-type galaxies offer a valu- 
able insight into their star formation and chemical evolution. 
We have used the VIMOS integral field unit to determine 
the gradients in spectroscopic features for 14 galaxies in the 
Shapley Supercluster. From these, we derived gradients in 
age, total metallicity and a-element abundance. We combine 
the Shapley gradients with revised values for our previous 



VIMO S observations (11 galaxies in Abell 3389: [Rawle et ah! 
l2008bl ). Our principal results are summarised below. 

• The average metallicity gradient is V[Z/H] = -0.13 ± 
0.04. 

• We find that the mean Vlog(age/Gyr) = -0.02 ± 0.06, 
but that ~40 per cent of galaxies have age gradients incon- 
sistent with zero (> 2 a). 

• Above a velocity dispersion of logcr ~ 2.1, galaxies ex- 
hibit a large scatter in metallicity gradient, suggesting that 
mergers with widely varying characteristics dominate their 
formation histories. 

• Galaxies with metal rich cores have steep negative 
V[Z/H] and have a preference for Vlog(age) > 0. Galax- 
ies with solar (or sub-solar) metallicity centres have V[Z/H] 
~ and Vlog(age) < 0. 

• There is a strong anti-correlation between the gradients 
of age and metallicity, with a slope of -0.69 ± 0.19. We 
have demonstrated that this relation cannot be attributed 
to correlated measurement errors. 

• Optical colours and colour gradients derived from the 
VIMOS spectroscopic age and metallicity are generally con- 
sistent with values obtained directly from imaging. Compar- 
ison between the spectroscopic and photometric B-R gradi- 
ents has a reduced ~ 3-4. 

A larger set of IFU observations, probing more of the 
galaxy parameter space (particularly towards lower mass), 
is required to explore the relationship between the gradi- 
ents further. Our results have successfully demonstrated the 
power of the integral field unit in measuring stellar popula- 
tion gradients for local cluster {z < 0.05) early-type galaxies. 
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Table Al. Median circular aperture corrections I(re) — I(re/3) 
for absorption line indices and kinematic/stellar population pa- 
rameters. Medians, calculated for the 25 galaxies with derived 
gradients, are in A, except CNi, CN2, Mgi and Mg2 which are 
quoted in mags. Ql and Q3 give the first and third quartiles of 
the data. 





median correction 


Ql 


Q3 


H(5A 


0.22 


± 


0.31 


-0.38 


0.74 


H(5F 


0.25 


± 


0.19 


-0.13 


0.50 


CNi 


-0.02 


± 


0.01 


-0.03 


0.01 


CN2 


-0.01 


± 


0.01 


-0.04 


0.00 


Ca4227 


-0.18 


± 


0.13 


-0.39 


0.04 


G4300 


-0.01 


± 


0.16 


-0.48 


0.28 


H7A 


0.22 


± 


0.17 


-0.36 


0.57 


H7F 


0.08 


± 


0.11 


-0.11 


0.44 


Fe4383 


0.13 


± 


0.19 


-0.46 


0.59 


Ca4455 


-0.11 


± 


0.06 


-0.25 


0.01 


Fe4531 


-0.07 


± 


0.11 


-0.32 


0.17 


Fe4668 


-0.81 


± 


0.20 


-1.57 


-0.36 


H/3 


0.07 


± 


0.08 


-0.06 


0.27 


Fe5015 


-0.12 


± 


0.18 


-0.40 


0.52 


Mgi 


-0.02 


± 


0.00 


-0.02 


-0.01 


Mg2 


-0.02 


± 


0.00 


-0.03 


-0.01 


Mgb5177 


-0.35 


± 


0.08 


-0.59 


0.01 


Fe5270 


-0.27 


± 


0.12 


-0.57 


-0.10 


Fe5335 


-0.34 


± 


0.12 


-0.49 


0.09 


Fe5406 


-0.11 


± 


0.07 


-0.25 


0.14 


Fe5709 


0.01 


± 


0.10 


-0.15 


0.19 


Fe5782 


-0.02 


± 


0.05 


-0.20 


0.09 


logo- 


-0.05 ± 0.01 


-0.06 


-0.01 


log{age/Gyr) 


-0.03 ± 0.05 


-0.18 


0.05 


[Z/H] 


-0.07 ± 0.04 


-0.19 


0.04 


[a/Fe] 


-0.02 ± 0.03 


-0.10 


0.02 
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APPENDIX A: APERTURE CORRECTION 

Comparison between spectroscopic absorption line strengths 
from different aperture sizes, requires a correction to ac- 
count for the gradient. Using the VIMOS IFU data, we have 
derived aperture corrections in the form I(re) — I(rc/3) for 
each index, and for log a, log(age), [Z/H] and [a/Fe]. Ta- 
ble lAll presents the median corrections for the 25 galaxies 
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with derived gradients (i.e. excluding D53 from A3389 and 
MGP0129, MGP1211 MGP2146, MGP2399, MGP4358 from 
Shapley) . 



APPENDIX B: CENTRAL VALUE AND 
GRADIENT DATA 

This paper presents VIMOS observations of 19 galaxies in 
the Shapley Supercluster core. Table IBll lists the central 
age, total metallicity and a-element over-abundace, as de- 
rived from spectra observed using an re/3 aperture. Table 
IB2I presents the stellar population gradients for 14 of the 
galaxies. Table lB3l lists two sets of B-R central colours and 
colour gradients for each galaxy: photometric from imaging; 
predicted from VIMOS stellar populations. 



APPENDIX C: REVISED RESULTS FOR 
RAWLE ET AL. (2008B) SAMPLE (A3389) 

The Shapley sample presented in this paper is combined 
with th e galaxies from Abell 3389 introduced in lRawle et al.l 
l|2008bl ). The raw data was re-analysed using the improved 
techniques described. Particularly important is the new and 
more stable implementation of the gradient computation, 
which requires 2-3 times fewer grid inversions and no lin- 
ear fitting to noisy stellar population data points. Table [CTI 
presents the central velocity dispersion and stellar popula- 
tion parameters (using re/3 apertures), while Table [C2l lists 
the revised gradients. 
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Table Bl. Central velocity dispersion and derived stellar populations for the 19 galaxies in the Shapley sample. Central parameters are 
calculated from an re/3 aperture. 



ID log a log(age/Gyr) [Z/H] [a/Fe] 



MGP0129 


2, 


.015 


± 


0, 


.026 


0, 


,867 


± 


0, 


,115 


-0, 


,024 


± 


0, 


,061 


0, 


,090 


± 


0, 


,089 


MGP1189 


2, 


.461 


± 


0, 


,007 


0, 


,945 


± 


0, 


,039 


0, 


,249 


± 


0, 


,030 


0, 


,240 


± 


0, 


,028 


MGP1195 


2, 


.063 


± 


0, 


,012 


0, 


,869 


± 


0, 


,068 


0, 


,204 




0, 


,042 


0, 


144 


± 


0, 


,038 


MGP1211 


1, 


.949 


± 


0, 


,050 


1, 


,013 


± 


0, 


,079 


0, 


,007 


± 


0, 


101 


-0. 


016 


± 


0, 


129 


MGP1230 


1 


.759 


± 


0, 


,029 


0, 


,734 


it 


0, 


,076 


-0, 


,171 


± 


0, 


,052 


0, 


133 


± 


0, 


,076 


MGP1440 


2, 


.273 


± 


0, 


,013 


1, 


,213 


± 


0, 


,049 


-0, 


,012 




0, 


,038 


0, 


,334 


± 


0, 


,037 


MGP1490 


1, 


.885 


± 


0, 


,028 


0, 


,639 


± 


0, 


,082 


0, 


,119 


± 


0, 


,082 


0, 


,186 


± 


0, 


,078 


MGP1600 


2, 


.473 


± 


0, 


,007 


1, 


,148 


± 


0. 


,026 


0, 


,274 


± 


0. 


,022 


0. 


,362 


± 


0. 


,018 


MGP1626 


2, 


.414 


± 


0, 


,012 


0, 


,491 


± 


0. 


111 


0, 


,660 


± 


0. 


,130 


0. 


,260 


± 


0. 


,042 


MGP1835 


2, 


.306 


± 


0, 


,009 


0, 


,874 




0. 


,048 


0, 


,475 


± 


0. 


,028 


0. 


,177 


± 


0. 


,020 


MGP1988 


2, 


.245 


± 


0, 


,014 


0, 


,778 


± 


0, 


,054 


0, 


,299 


± 


0, 


,036 


0, 


,303 


± 


0, 


,036 


MGP2083 


2, 


.387 


± 


0, 


,006 


1, 


,014 


± 


0, 


,027 


0, 


,160 


± 


0, 


,028 


0, 


,222 


± 


0, 


,021 


MGP2146 


2, 


.084 


± 


0, 


,015 


0, 


,836 


± 


0, 


,080 


0, 


,109 


± 


0, 


,044 


0, 


,225 


± 


0, 


,045 


MGP2399 


1, 


.906 


± 


0, 


,042 


0, 


,662 


± 


0, 


,098 


0, 


,073 


± 


0, 


106 


0, 


,051 


± 


0, 


,088 


MGP2437 


2, 


.271 


± 


0, 


,014 


1, 


,412 


± 


0, 


,053 


0, 


,100 




0, 


,043 


0, 


,382 


± 


0, 


,035 


MGP2440 


2, 


.356 


± 


0, 


,007 


1, 


,043 




0, 


,019 


0, 


,096 


± 


0, 


,022 


0, 


,206 


± 


0, 


,020 


MGP3971 


2, 


.261 


± 


0. 


,022 


1, 


,144 


± 


0. 


,051 


0, 


,117 


± 


0. 


,051 


0. 


,349 


± 


0. 


,044 


MGP3976 


2, 


.402 


± 


0. 


,009 


1, 


,118 


± 


0. 


,031 


0, 


,194 


± 


0. 


,030 


0. 


,358 


± 


0. 


,025 


MGP4358 


1, 


.832 


± 


0. 


,070 


0, 


,659 


± 


0. 


,153 


0, 


,237 


± 


0. 


,148 


-0. 


,188 


± 


0. 


,120 



Table B2. Velocity dispersion and stellar population gradients for 14 galaxies in the Shapley sample. 



ID Vlogo- Vlog(age/Gyr) V[Z/H] V[q/Fc] 



MGP1189 


-0, 


,137 


± 


0, 


,040 





.018 


± 


0. 


,119 


-0, 


,133 


± 


0, 


,084 


-0, 


,067 


± 


0, 


,102 


MGP1195 


-0, 


,094 


± 


0, 


,057 





.013 


± 


0, 


,201 


-0, 


,294 


± 


0, 


,141 


-0, 


,066 


± 


0, 


,123 


MGP1230 


-0, 


,073 


± 


0, 


,130 


-0, 


.413 


± 


0, 


,195 


0, 


,066 


± 


0, 


,185 


-0, 


,310 


± 


0, 


,167 


MGP1440 


-0, 


152 


± 


0, 


,058 


-0, 


,259 


± 


0, 


,221 


-0, 


136 


± 


0, 


,144 


-0, 


,038 


± 


0, 


,148 


MGP1490 


-0, 


,038 


± 


0, 


,031 


-0, 


,195 


± 


0, 


,250 


0, 


,050 


± 


0, 


,201 


-0, 


,245 


± 


0, 


,244 


MGP1600 


-0, 


,071 


± 


0, 


,028 


0, 


.022 


± 


0, 


,114 


-0, 


,166 


± 


0, 


,095 


0, 


102 


± 


0, 


,075 


MGP1626 


-0. 


,145 


± 


0, 


,046 





.484 


± 


0. 


,120 


-0. 


,508 


± 


0, 


,126 


-0. 


,042 


± 


0, 


,074 


MGP1835 


-0. 


,088 


± 


0, 


,026 


-0, 


.029 


± 


0. 


,089 


-0. 


,172 


± 


0, 


,054 


0. 


,055 


± 


0, 


,052 


MGP1988 


-0. 


,179 


± 


0, 


,039 





.031 


± 


0. 


,209 


-0. 


,148 


± 


0, 


,128 


-0. 


,090 


± 


0, 


,090 


MGP2083 


-0. 


,046 


± 


0, 


,029 


-0, 


.251 


± 


0. 


,088 


-0. 


,093 


± 


0, 


,085 


-0. 


,045 


± 


0, 


,064 


MGP2437 


-0. 


,232 


± 


0, 


,112 


-0, 


.523 


± 


0. 


,153 


0. 


,201 


± 


0, 


,110 


-0. 


,014 


± 


0, 


,120 


MGP2440 


-0. 


,036 




0, 


Oil 


-0, 


.265 


± 


0, 


,085 


0, 


,025 


± 


0, 


,083 


-0. 


,058 


± 


0, 


,064 


MGP3971 


-0, 


,088 


± 


0, 


,010 


0, 


.306 




0, 


,201 


-0, 


,310 




0, 


,137 


-0, 


,034 


± 


0, 


,121 


MGP3976 


-0, 


,198 


± 


0, 


,038 


-0, 


.211 


± 


0, 


,174 


0, 


,147 


± 


0, 


,138 


-0, 


,014 


± 


0, 


,113 
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Table B3. Central B-R colour {rc/3 aperture) and colour gradient for each galaxy in the Shapley sample. Photometric measurements 
are from WFI B- and _R-band imaging, in Vega magnitudes and corrected for galactic extinction. Colour gradients are calculated using 
the same annuli as the spectroscopic gra dients. Predicted values are derived from the VIMOS spectroscopic age and metallicity, via the 
evolutionary stellar population models of iMarastorJ ||2005|) 



ID Photometric VIMOS predicted 

central B-R VB-R central B-R VB-R 



MGP0129 


1, 


.62 






- 






1, 


,51 












MGP1189 


1 


.63 


-0. 


.15 


± 


0. 


,04 


1 


,62 


-0. 


,04 


± 


0. 


,04 


MGP1195 


1, 


.57 


-0. 


.08 


± 


0. 


,02 


1, 


,59 


-0. 


,10 


± 


0. 


,01 


MGP1211 


1, 


.61 












1, 


,56 












MGP1230 


1, 


.55 


-0, 


.05 


± 


0. 


,03 


1, 


,43 


-0. 


,05 


± 


0. 


,03 


MGP1440 


1, 


.64 


-0. 


.19 


± 


0. 


,03 


1, 


,65 


-0. 


18 


± 


0. 


,02 


MGP1490 


1, 


.54 


-0. 


.14 


± 


0. 


,01 


1, 


,50 


-0. 


,04 


± 


0. 


,04 


MGP1600 


1, 


.69 


-0. 


.05 


± 


0. 


,02 


1, 


,72 


-0. 


,04 


± 


0. 


,01 


MGP1626 


1, 


.70 


-0. 


.11 


± 


0. 


,02 


1, 


,64 


-0. 


,02 


± 


0. 


,03 


MGP1835 


1 


.65 


-0. 


.05 


± 


0. 


,03 


1, 


,68 


-0. 


,06 


± 


0. 


,02 


MGP1988 


1 


.71 


-0. 


.14 


± 


0. 


,01 


1, 


,61 


-0. 


,05 


± 


0. 


,03 


MGP2083 


1 


.75 


-0. 


.20 


± 


0. 


,02 


1, 


,62 


-0. 


12 


± 


0. 


,01 


MGP2146 


1, 


.52 












1, 


,55 












MGP2399 


1, 


.52 












1, 


,49 












MGP2437 


1 


.55 


-0. 


.07 


± 


0. 


,01 


1, 


,79 


-0. 


19 


± 


0. 


,07 


MGP2440 


1 


.58 


-0. 


.07 


± 


0. 


,02 


1, 


,61 


-0. 


,07 


± 


0. 


,01 


MGP3971 


1 


.74 


-0. 


.11 


± 


0. 


,07 


1, 


,66 


+0. 


,03 


± 


0. 


,06 


MGP3976 


1, 


.63 


-0. 


.10 


± 


0. 


,02 


1, 


,68 


-0. 


,02 


± 


0. 


,04 


MGP4358 


1 


.59 












1, 


,55 













Table CI. Revise d central velocity dispersion and stellar populations for the 12 galaxies in the A3389 sample (observations from 
iRawle et al.ll2008bl) . Central parameters are calculated from an rc/3 aperture. 



ID log a log(age/Gyr) [Z/H] [a/Fc] 



D40 


2. 


.273 


± 


0. 


,012 


1, 


.006 


± 


0. 


.038 


0. 


.230 


± 


0. 


,042 


0. 


,369 


± 


0. 


.039 


D41 


2. 


.315 


± 


0. 


,005 


1, 


.018 


± 


0. 


.021 


0. 


.347 


± 


0. 


,018 


0. 


137 


± 


0. 


.013 


D42 


2. 


.455 


± 


0. 


,005 


1, 


.008 


± 


0. 


.019 


0. 


.339 


± 


0. 


,017 


0. 


,324 


± 


0. 


.012 


D43 


2. 


.194 


± 


0. 


,009 


1, 


.044 


± 


0. 


.037 


0. 


.233 


± 


0. 


,037 


0. 


,073 


± 


0. 


.030 


D44 


2. 


,406 


± 


0. 


,009 


0, 


.928 


± 


0. 


.047 


0. 


.380 


± 


0. 


,025 


0. 


,239 


± 


0. 


.024 


D48 


2. 


,216 


± 


0. 


,009 


1, 


.077 


± 


0. 


.030 


0. 


.241 


± 


0. 


,030 


0. 


,083 


± 


0. 


.025 


D49 


2. 


,306 


± 


0. 


Oil 


1, 


.331 


± 


0. 


.042 


0. 


.101 


± 


0. 


,035 


0. 


,391 


± 


0. 


.026 


D52 


2. 


,446 


± 


0. 


,005 


1, 


.310 


± 


0. 


.019 


0. 


.128 


± 


0. 


,016 


0. 


,292 


± 


0. 


.012 


D53 


2. 


,077 


± 


0. 


,012 


0, 


.642 


± 


0. 


.072 


0. 


.354 


± 


0. 


,048 


0. 


182 


± 


0. 


.042 


D60 


2. 


,317 


± 


0. 


,003 


0, 


.805 


± 


0. 


.030 


0. 


.308 


± 


0. 


,015 


0. 


159 


± 


0. 


.013 


D61 


2. 


,168 


± 


0. 


,011 


0, 


.652 


± 


0. 


.064 


0. 


.340 


± 


0. 


,045 


0. 


,076 


± 


0. 


.037 


D89 


2. 


,416 


± 


0. 


,006 


1, 


.033 


± 


0. 


.022 


0. 


.304 


± 


0. 


,020 


0. 


,217 


± 


0. 


.015 



Table C2. Revised velocity dispersion and stellar population gradients for 11 galaxies in the A3389 sample (D53 is excluded). 



ID Vlogo- Vlog(age/Gyr) V[Z/H] V[a/Fe] 



D40 


0.075 


± 


0. 


.085 


-0.106 


± 


0.184 


-0.329 


± 


0.115 


0.258 


± 


0. 


.128 


D41 


0.037 


± 


0. 


.035 


0.666 


± 


0.062 


-0.407 


± 


0.045 


0.113 


± 


0. 


.034 


D42 


-0.076 


± 


0. 


.008 


-0.138 


± 


0.103 


-0.006 


± 


0.056 


-0.058 


± 


0. 


.041 


D43 


-0.072 


± 


0. 


.023 


0.297 


± 


0.181 


-0.266 


± 


0.131 


0.012 


± 


0. 


.098 


D44 


-0.097 


± 


0. 


.018 


0.467 


± 


0.136 


-0.224 


± 


0.105 


-0.254 


± 


0. 


.083 


D48 


-0.090 


± 


0. 


.049 


-0.164 


± 


0.105 


0.156 


± 


0.094 


-0.262 


± 


0. 


.076 


D49 


-0.036 


± 


0. 


.006 


-0.047 


± 


0.188 


-0.053 


± 


0.135 


-0.252 


± 


0. 


.121 


D52 


0.002 


± 


0. 


.012 


-0.159 


± 


0.048 


-0.011 


± 


0.043 


-0.195 


± 


0. 


.028 


D60 


-0.069 


± 


0. 


.011 


-0.168 


± 


0.066 


-0.298 


± 


0.053 


-0.177 


± 


0. 


.047 


D61 


-0.046 


± 


0. 


.064 


-0.026 


± 


0.184 


-0.054 


± 


0.144 


-0.238 


± 


0. 


.112 


D89 


-0.268 


± 


0. 


.042 


0.139 


± 


0.070 


-0.410 


± 


0.064 


-0.685 


± 


0. 


.045 



